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Using 88.9 million BB events collected by the BABAR detector at the T(4S), we measure the 
branching fraction for the radiative penguin process B — » X a "/ from the sum of 38 exclusive fi- 
nal states. The inclusive branching fraction above a minimum photon energy _E 7 >1.9 GeV is 
B{b -> S7) = {3.27±OA8(stat.)to.ll(syst.)1:o^i(theory)) x 10~ 4 . We also measure the isospin asym- 
metry between B~ -> X 3 al and B° -> X sS -y to be A _ = -0.006 ± 0.058(stat.) ± 0.009(syst.) ± 
0.024(5° /B~). The photon energy spectrum is measured in the B rest frame, from which moments 
are derived for different values of the minimum photon energy. We present fits to the photon spec- 
trum and moments which give the heavy-quark parameters mj and fi^. The fitted parameters are 
consistent with those obtained from semileptonic B — > X c lv decays, and are useful inputs for the 
extraction of \ V u b\ from measurements of semileptonic B — > X u lv decays. 

PACS numbers: 13.35.Dx, 14.60.Fg, 11.30.Hv 



I. INTRODUCTION 

Radiative decays involving the flavor-changing neutral 
current transition b — > s are described in the Standard 
Model primarily by a one-loop radiative penguin dia- 
gram containing a top quark and a W boson. Calcula- 
tions of this rate in the Standard Model have now been 
completed to next-to- leading order in a s (NLO), with 
a predicted branching fraction B(b — > sj) = (3.57 ± 



0.30) x 10~ 4 for E 7 > 1.6 GeV 



QBE), 



which is 



consistent with the current experimental world average 
B(b -> sj) = (3.52±g;lg) x 10" 4 0. Calculations of 
next-to-next-to- leading order (NNLO) corrections are in 
progress [Ej]. Additional contributions to the loop from 
new physics, e.g. a charged Higgs boson, could change 
the s 7 rate [1 0, @, ll D3 • 

The photon energy spectrum in B — > X s "f provides 
access to the distribution function of the b quark inside 
the B meson ^l). The knowledge of this shape function 
is a crucial input in the extraction of \V u h\ from inclusive 
semileptonic B — » X u tv measurements 0,0,0,0,0, 
Il7j . We fit the spectrum to two recent predictions, one 
using a combination of the operator product expansion 
(OPE) coupled to soft collinear effective theory 0, 0, 
0, llTl ITsl llflj in the shape function mass scheme, and 
the other using a full OPE approach in the kinetic mass 
scheme po| . 

We also present measurements in the B rest frame of 
the first, second and third moments of the photon energy 
spectrum for five different minimum energies, E 1 > 1.90, 
2.00, 2.09, 2.18 and 2.26 GeV. The heavy quark param- 
eters mi, and /n^, which describe the effective b quark 
mass and kinetic energy inside the B meson, can be 
determined either from fits to these moments [20j, or 
from the fits to the spectrum. We compare the fitted 
parameters with those obtained from the lepton energy 



and hadronic mass moments measured in semileptonic 
B -> X c lv decays I^tI . Previous measurements of the in- 
clusive branching fraction have used two different meth- 
ods. In the fully inclusive method |22l I23L |24| the photon 
energy spectrum is measured without reconstructing the 
X s system, and backgrounds are suppressed using in- 
formation from the rest of the event. When measuring 
the E 1 spectrum inclusively at the T(4S) the shape of 
the spectrum has a large contribution from the 50 MeV 
calorimeter energy resolution, and from the motion of the 
B meson in the T(4S) rest frame. 

The semi-inclusive method [2^, 0, |2(| uses a sum of 
exclusive final states in which possible X s systems are 
combined with the photon, and kinematic constraints of 
T(4S) production are used to suppress backgrounds. We 
have chosen this method and made several improvements 
over the previous analyses. The number of final states of 
the X s system has been increased to 38 by the inclusion 
of states with two 7r°s, 77 mesons, and three kaons, and 
the X s mass range has been increased to 0.6 — 2.8 GeV. 
Candidates with correctly reconstructed X s systems are 
treated as signal, whereas incorrectly reconstructed sys- 
tems, referred to as "cross-feed", are treated as back- 
ground. This method allows us to make a measurement 
of the branching fraction as a function of the hadronic 
mass, M(X S ). The M(X S ) spectrum is converted into a 
high resolution photon energy spectrum using the kine- 
matic relationship for the decay of a B meson of mass 
M B : 



Ej = 



M% - M{X s f 
2Mb 



(1) 



where E 7 is the photon energy in the B rest frame which 
has a resolution of 1-5 MeV. 



II. DETECTOR AND DATA 



*Also with Universita di Perugia, Dipartimento di Fisica, Perugia, 
Italy 

tAlso with Universita della Basilicata, Potenza, Italy 
t Deceased 



The results presented in this paper are based on data 
collected with the BABAR detector :27] at the PEP-II 
asymmetric-energy e + e~ collider located at the Stanford 
Linear Accelerator Center. The data sample has an in- 
tegrated luminosity of 81.9 fb _1 , corresponding to 88.9 
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million BB pairs recorded at the T(AS) resonance ("on- 
peak", center-of-mass energy y/s — 10.58 GeV). An ad- 
ditional 9.6 fb _1 were recorded about 40 MeV below this 
energy ("off-peak"), for the study of continuum back- 
grounds in which a light or charm quark pair is produced. 

The asymmetric beam configuration in the laboratory 
frame provides a boost of fi*f = 0.56 to the T(45). This 
results in the high energy photons from b — > sj de- 
cays having energies between 1.5 and 4.5 GeV in the 
laboratory frame. Photons are detected and their ener- 
gies measured by a CsI(Tl) electromagnetic calorimeter 
(EMC). The energy scale of the calorimeter crystals is 
determined by radioactive source and Bhabha scattering 
calibrations, and the energy leakage of photon showers 
is corrected using ir° decays. The photon energy resolu- 
tion is measured with symmetric ir° and i] decays to be 
a E /E = {2.3/ E( GeV) 1 / 4 © 1.4} %, where the terms are 
added in quadrature. The measured 7r° mass resolution 
is between 6 and 7 MeV for momenta below 1 GeV in 
the laboratory frame. 

Charged particles are detected and their momenta 
measured by the combination of a silicon vertex tracker, 
consisting of five layers of double-sided detectors, and 
a 40-layer central drift chamber, both operating in the 
1.5 T magnetic field of a solenoid. The transverse mo- 
mentum resolution for the tracking system is (J Pt /pt — 
0.0013pT ffi 0.0045, where pt is measured in GeV. 

Charged particle identification is provided by the av- 
erage energy loss (cLE / dx) in the tracking devices and by 
an internally reflecting ring-imaging Cherenkov detector 
(DIRC). The dE/dx resolution from the drift chamber 
is typically 7.5% for pions. The Cherenkov angle res- 
olution of the DIRC is measured to be 2.4mrad, which 
provides more than 3<r separation between charged kaons 
and pions up to a momentum of 3 GeV. 

III. X s SIGNAL MODEL AND BACKGROUNDS 

The B — > X s j signal includes resonant and non- 
resonant X s states, but S-wave states are forbidden by 
angular momentum conservation. It is known experimen- 
tally that the mass region M(X S ) < 1.1 GeV is domi- 
nated by the K* (892) resonance 28] . In the higher-mass 
region there is evidence for the If i (1270) and if|(1430) 
resonances [2^, but these only account for about 16 % of 
the inclusive rate. Theoretical predictions for exclusive 
decays to higher K* resonances also account for less than 
half of the inclusive rate [3(j ■ 

The sum of many broad resonances can be modeled 
by an inclusive distribution. In designing our analysis 
we have used an inclusive calculation from Kagan and 
Neubert [j| . We follow their prescription and replace the 
inclusive model in the region M(X S ) < 1.1 GeV with an 
equivalent amount of exclusive B — > if* (892)7. The Ka- 
gan and Neubert calculation has two empirical parame- 
ters, mb and Ai, which are related to the mean and width 
of the spectrum (they are similar to the OPE parameters 
TOfc and — /i 2 ). The default parameters were originally 



chosen to be mj, = 4.65 GeV and Ai = —0.52 GeV 2 , but 
eventually we fit for these parameters using our measured 
spectrum. In the inclusive region M(X S ) = 1.1 — 2.8 GeV 
the fragmentation of the X s system into hadrons is sim- 
ulated using JETSET [3l| . The response of the detector 
is modeled using GEANT4 

Most of the background in this analysis arises from 
continuum production of a high energy photon, either by 
initial state radiation, or from the decays of ir° and r\ 
mesons produced in light-quark jets. We combine event 
shape information into a neural network and use the out- 
put to remove most of this background. The ir° and r\ 
contributions are further suppressed by vetoes on combi- 
nations of the high energy photon with another photon 
in the event which have a mass consistent with a 7r° or 
V- 

Backgrounds from hadronic 6 — > c decays are impor- 
tant for M(X S ) > 1.8 GeV. There are two contributions 
to this background: combinatorial final states containing 
particles from both B decays, and incomplete final states 
where the particles all come from the decay of one B, but 
one or more of the decay products is missing. In the case 
when only one low energy photon is missing from a B 
decay to a final state containing a n°, e.g. B — > D^*'p~, 
the events tend to peak in the signal region, since only 
a small amount of energy is missing from the final state. 
At high hadronic masses this background becomes com- 
parable to the expected signal yield. 

Backgrounds from charmless hadronic B decays give 
a small contribution over the whole X s mass range, but 
they include a poorly understood component from B — > 
X s tt° decays which can peak in the signal region. These 
backgrounds are modeled as a sum of the measured and 
yet unmeasured charmless decay modes. There is also a 
cross-feed background from mis-reconstructed B — > X s j 
decays which is discussed in detail in Section [3 The 
requirement of a positively identified kaon removes B — > 
X d j decays. 

IV. B MESON RECONSTRUCTION 

We reconstruct the X s states in 38 decay modes and 
their charge conjugates. This includes 22 final states with 
a kaon and one to four pions, where at most two of the 
pions are 7r°s, 10 states with a kaon, an r\ and up to two 
pions, and 6 states with three kaons plus at most one 
pion. The kaons can be either K~ or K®. A full list 
of the modes can be found in Table |nj According to 
our signal model these modes represent 55 % of the total 
inclusive rate in the region M(X S ) = 1.1 — 2.8 GeV. The 
X s modes that we do not reconstruct are referred to as 
"missing fractions" . 

Neutral kaons are reconstructed as K® — > tt + tt^ candi- 
dates with an invariant mass within 9 MeV of the nomi- 
nal K° mass j3^| , and a transverse flight distance greater 
than 2 mm from the primary event vertex. We do not re- 
construct K ° — * 7r°7r° because of its low efficiency, and we 
do not reconstruct K® because we cannot directly mea- 
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sure its energy. Charged kaons are identified using in- 
formation from the DIRC and the tracking devices. The 
remaining tracks are considered to be from charged pi- 
ons. Both charged and neutral kaons are required to have 
momenta > 0.7 GeV in the laboratory frame. Above this 
threshold the rate for charged pions to be mis-identified 
as kaons is < 2%. 

Neutral pions are reconstructed from pairs of photons, 
each with an energy > 30 MeV. For 7r° candidates in 
the mass interval 117 and 150 MeV, a fit is performed to 
improve the momentum resolution. To reject combinato- 
rial background, charged and neutral pions are required 
to have momenta > 0.5, 0.3 or 0.2 GeV in the laboratory 
frame for states with 1, 2 or > 3 pions, respectively. 

The r) mesons are reconstructed from pairs of photons 
with energies > 50 MeV. For r\ candidates in the mass 
interval 520 and 580 MeV, a fit is performed to improve 
the momentum resolution. The rj mesons are required to 
have momenta > 0.7 GeV in the laboratory frame. We 
do not explicitly reconstruct the modes rj 7r + 7r~7r° and 
r\ — * 7r 7r°7r°, but the former decays are included in the 
final states with a kaon and > 3 pions. 

The reconstructed X s system is combined with a high- 
energy photon to form a B meson. The photon is de- 
tected as an isolated energy cluster in the calorimeter, 
with shape consistent with a single photon, and energy 
E* > 1.8 GeV in the e + e~ center-of-mass (CM) frame. 
A veto is applied to high energy photons that, combined 
with another photon, form either a 7r° within the mass 
range 117-150 MeV or an r\ within the mass range 524- 
566 MeV. In the higher-mass region M(X S ) > 2.0 GeV, 
we improve the rejection of B — ► D^p~ background by 
opening up the ir° mass window to 106-162 MeV and 
applying a veto if the ir° forms a p~ meson in the mass 
range 400-1200 MeV when combined with a charged pion. 

We remove 85 % of the continuum background by se- 
lections on the angle, 0^, between the thrust axis of the 
B meson candidate and the thrust axis of all the other 
particles of the event, requiring | cos#Jn| < 0.80, and the 
angle, 9 B , between the B candidate and the beam axis, 
requiring \cos9* B \ < 0.80, both defined in the e + e~~ CM 
system. We then use a neural network to combine infor- 
mation from a set of event-shape variables. The inputs 
to the neural net include 9, 



tween the measured and expected energies of the B can- 



didate, AE = E% - (y/s/2), 



where E* B and p* B are the 



T' V B> 



R2, the ratio of the sec- 
ond to zeroth Fox- Wolfram moments [34j. and R2' and 
9' T , which are defined in the primed frame obtained by 
removing the high energy photon and boosting the rest 
of the event into its CM frame. These last two variables 
discriminate against background from initial state radi- 
ation. The inputs to the neural net also include a set 
of 18 energy flow cones each covering an angle of 10° 
about the reconstructed B direction in the e + e~ CM 
system. The neural network selection is tightened above 
M(X S ) = 1.1 GeV, and again above M(X S ) = 2.0 GeV 
and M(X S ) = 2.4 GeV, to take account of the increasing 
background as a function of hadronic mass. 

The identification of B — > V s 7 decays makes use of 
two kinematic variables: the beam-energy substituted 

and the difference be- 



energy and momentum of the B candidate in the CM 
frame, and y/s is the total CM energy. When calculating 
wies, the value of p* B is corrected for the tail of the high 
energy photon response function of the EMC by scaling 
the measured E* to the value that gives AE = 0, the 
value expected for true signal. 

Within an initial selection toes > 5.22 GeV and 
\AE\ < 0.40 GeV we reconstruct two candidates per 
event on average. In events where more than one B 
candidate is reconstructed we select the best candidate 
using the smallest value of |Ai?|. This technique is > 
90% efficient when the true B — > X s j decay is among 
the reconstructed candidates. The \AE\ distribution 
has a resolution of about 0.05 GeV with a radiative tail 
on the low side. For the best candidate we require 
\AE\ < 0.10 GeV for M(X S ) < 2.0 GeV, and tighten 
this to \ AE\ < 0.08 GeV for M(X S ) = 2.0-2.4 GeV and 
\AE\ < 0.07 GeV for M(X S ) = 2.4 - 2.8 GeV. These 
selections are optimized to give the best statistical sig- 
nificance for the signal in each X s region. We then fit the 
rn.ES distribution between 5.22 and 5.29 GeV to extract 
the signal yield. 




1.0 1.4 1.8 2.2 

M (X s ) (GeV) 



FIG. 1: Efficiency for correctly reconstructing a signal event 
in one of the 38 final states as a function of hadronic mass. 
Note that this efficiency does not include the missing fractions 
of B — > X S J final states. 



V. SIGNAL EFFICIENCY AND CROSS-FEED 

The signal efficiency is determined from generated 
Monte Carlo events which are produced and correctly 
reconstructed in one of the 38 final states. It does not 
include the missing B — ► V s 7 final states which are dis- 
cussed in Section IV11I We generate equal numbers of 
B~ and B° decays and assume isospin symmetry. The 
production of K~ and K° is equal, and the branching 



fraction for the K 1 



7r + 7r decay is included in the 



mass, m ES = (y r s/2) 2 



Pb 



Monte Carlo generator. For the isospin asymmetry mea- 
surement we note that the efficiency for reconstructing 
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B~ decays is lower than for B° decays by almost a fac- 
tor of two. This difference results from a combination 
of the different distributions of X s final states and the 
different efficiencies for K® and K~ , 7r~ and ir°. 

The efficiency is a rapidly varying function of hadronic 
mass. In the K* region, M(X S ) < 1.1 GcV, the efficiency 
is about 15%, dominated by the high efficiency for re- 
constructing the K — » K~n + mode. For M(X S ) > 
1.1 GeV the efficiency decreases from 5 % to 1.5 % as the 
hadronic mass increases, as shown in Figure ^ There 
are two main reasons for the mass dependence: the mul- 
tiplicity of the final state particles increases with mass, 
and the angular correlation between them decreases. In 
addition, there are steps in efficiency at M(X S ) = 1.1, 2.0 
and 2.4 GeV because the selection criteria are tightened 
as the levels of background increase with mass. 

The efficiencies that we obtain from the signal samples 
are corrected for small differences in detection efficiencies 
between data and Monte Carlo events which are deter- 
mined using control samples. The tracking efficiency is 
reduced by (0.8 ± 2.0)% per track, and the K® efficiency 
by (2 ± 3)%. The photon efficiency is not adjusted, but 
is assigned a 2.5% error per photon. The effect of these 
adjustments is to reduce the overall signal efficiency by 
(1.9 ± 6.0)%. The final state distributions are also ad- 
justed by re-weighting the signal Monte Carlo events to 
match the distributions in data (see Section IVIIfl . 

Monte Carlo events that are reconstructed in a dif- 
ferent final state, or with the wrong hadronic mass, are 
treated as a cross-feed background. The main sources of 
cross-feed are events from missing final states which are 
reconstructed in one of the 38 final states, and events 
from the 38 final states in which one of the final state 
particles from the X s is undetected and replaced by a 
low momentum particle from the other B. There is a 
small contribution from events with multiple candidates 
in which the true candidate is rejected in favour of an- 
other candidate with a smaller value of |Ai?|. 

The cross-feed background is fitted together with the 
hadronic B decay backgrounds, and allowed to have a 
component that peaks in the signal region. The amount 
of cross- feed increases as a function of hadronic mass, and 
is largest for the high multiplicity final states. As part 
of our systematic studies we vary the definition of cross- 
feed, transferring events between the signal and cross- 
feed samples. This gives changes in background and sig- 
nal efficiencies which compensate each other within 1 %. 



VI. FITTING 

To extract the signal yield we fit the toes distribu- 
tion of the data with a combination of a Crystal Ball 
function [35| for signal, a Novosibirsk function (33| for 
the peaking backgrounds, and ARGUS functions |3j| for 
the combinatorial backgrounds. The final data fits use 
Toes shapes derived from fits to three sets of Monte Carlo 
samples: signal, continuum background, and BB back- 
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FIG. 2: Fits to the M{X S ) = 1.4 - 1.5 GeV bin for (a) 
correctly reconstructed signal Monte Carlo events; (b) sim- 
ulated continuum background; (c) simulated cross-feed and 
hadronic B decay backgrounds, where the contributions from 
the peaking (dotted) and combinatorial (dashed) backgrounds 
are shown separately; (d) on-peak data, where the contri- 
butions from the signal (dotted-dashed) , continuum back- 
ground(dotted), peaking B background (long-dashed), and 
combinatorial B background (dashed) are shown separately. 



ground which includes cross-feed. In all cases we use an 
unbinned maximum likelihood fit. 

Figure |2Ja) shows the fit to the toes distribution of 
correctly reconstructed signal Monte Carlo events with 
M(X S ) — 1.4 — 1.5 GeV. We find no significant varia- 
tions in the parameters of the Crystal Ball shape over the 
M(X S ) range, so the values used in the fits to the data 
are fixed to the weighted average of the results over the 
full M(X S ) range: a width a = (2.81 ± 0.05) MeV, a tail 
parameter a = 2.17 ± 0.12, and a slope n = 0.99 ± 0.19. 

Figure |5{b) shows a fit of an ARGUS function to con- 
tinuum Monte Carlo events in the same bin in M(X S ). 
The amount of continuum background increases with 
M(X S ), and the shape parameter of the ARGUS func- 
tion is a rapidly varying function of M(X S ). To cross- 
check our understanding of the continuum background 
we also fit the off-peak data sample. While the overall 
variations are well reproduced, we find small systematic 
differences in both the shape parameter and the normal- 
ization which can only be partially accounted for by the 
difference between the on- and off-peak center-of-mass 
energies. We fix the continuum ARGUS shapes and nor- 
malizations in the fit to the data to the values from the 
continuum Monte Carlo, which vary from bin to bin, and 
are adjusted for small offsets observed in comparison to 
off-peak data. The difference in the normalization be- 
tween off-peak data and continuum Monte Carlo is 2.8 %. 
These offsets are considered to be part of the systematic 
errors. 

Figure |2Ic) shows a fit to the sum of the Monte Carlo 
predictions for the cross-feed and the hadronic B decay 
backgrounds using an ARGUS function plus a peaking 
Novosibirsk function. In the region M(X S ) < 1.8 GeV 
the largest contribution comes from the cross-feed, with 
only a small contribution from charmless hadronic B de- 



8 



cays. In the region M(X S ) > 1.8 GeV the hadronic b — > c 
background increases rapidly. The peaking background 
shape is broader than the signal shape, reflecting the less- 
peaked behavior of the backgrounds. The shape of the 
Novosibirsk function is determined to have a — 5.0 MeV 
and t = —0.295, from a fit to the simulated BB back- 
grounds over the full M(X S ) range. 

Figure E^d) shows the fit to the on-peak data in the 
same mass bin. In this fit the signal yield and the shape 
and normalization of the combinatorial BB background 
function are allowed to vary. The signal shape is fixed, 
and the continuum ARGUS function and peaking BB 
background shapes are fixed. The peaking background 
yield is obtained from the fit shown in Figure Gfc). Ta- 
ble Ogives the fitted signal yields, the peaking BB back- 
ground yields, and, as a measure of the goodness of the 
maximum likelihood fit, the \ 2 P er degree of freedom. 

Figure 13 shows the on-peak data fit to the full M(X S ) 
range, which gives the yield in the last row of Tabled In 
this fit the continuum ARGUS shape, and the peaking 
BB background are taken from fits to the full M(X S ) 
range of the simulated continuum and BB samples. 

The fit procedure has been validated with Monte Carlo 
studies to check for biases in the fitting method. System- 
atic errors from the fitting method, including variations 
in the fixed parameters in the fits, are discussed in Sec- 
tion Eini 

TABLE I: Signal yields from the fits to the on-peak data and 
the xVdof from the fits. Also given are the peaking back- 
ground yields from fits to cross-feed and BB Monte Carlo. 



M(X a ) Data Signal Data Fit Peaking Bkg 
(GeV) yield (events) xVdof yield (events) 



0.6-0.7 


6.5 


± 


7.7 


2.2 


0.9 


± 


3.3 


0.7-0.8 


5.6 


± 


14.1 


0.8 


2.7 


± 


6.4 


0.8-0.9 


416.2 


± 


23.2 


1.5 


24.2 


± 


8.5 


0.9-1.0 


355.6 


± 


24.9 


0.9 


22.7 


± 


10.8 


1.0-1.1 


51.3 


± 


19.0 


1.0 


14.4 


± 


13.7 


1.1-1.2 


33.2 


± 


12.9 


1.2 


7.4 


± 


6.0 


1.2-1.3 


83.2 


± 


15.7 


1.1 


9.4 


± 


7.5 


1.3-1.4 


101.5 


± 


16.8 


0.8 


0.8 


± 


8.6 


1.4-1.5 


72.0 


± 


15.8 


0.8 


15.3 


± 


9.1 


1.5-1.6 


82.4 


± 


16.5 


1.1 


16.1 


± 


11.3 


1.6-1.7 


66.1 


± 


16.9 


1.0 


5.3 


± 


11.7 


1.7-1.8 


54.6 


± 


16.5 


1.3 


5.6 


± 


13.1 


1.8-1.9 


76.6 


± 


18.2 


1.1 


19.1 


± 


13.7 


1.9-2.0 


13.5 


± 


19.5 


1.1 


21.3 


± 


14.1 


2.0-2.2 


47.5 


± 


21.8 


0.7 


19.4 


± 


16.8 


2.2-2.4 


52.1 


± 


24.0 


0.7 


39.5 


± 


21.0 


2.4-2.6 


44.7 


± 


25.6 


0.8 


46.8 


± 


20.8 


2.6-2.8 


-6.2 


± 


31.9 


1.0 


81.0 


± 


26.1 


0.6-2.8 


1513.0 


± 


85.1 


1.2 


464.2 


± 


68.3 




m ES (GeV) 

FIG. 3: On-peak data fit to the full M(X S ) range, M(X S ) = 
0.6 — 2.8 GeV bin, where the contributions from the signal 
shape (dotted-dashed), fixed continuum ARGUS shape (dot- 
ted), peaking background shape (long-dashed), and combina- 
torial BB background shape (dashed) are shown separately. 



VII. X s FRAGMENTATION AND MISSING 
FRACTIONS 

The fragmentation of the X s system into hadronic final 
states has been modeled using JETSET .3jJ. We check 
this fragmentation by comparing the observed data yields 
in the range M(X S ) = 1.1 - 2.8 GeV with the yields 
expected from the signal Monte Carlo. We do a detailed 
study by splitting the data and Monte Carlo samples 
into 10 different categories, each containing two to ten of 
our selected final states. The measured ratios of fitted 
data signal events to reconstructed Monte Carlo signal 
events in each category are given in Table [H] with their 
statistical errors. 

We note that the rates for the Kit modes are smaller 
than the JETSET prediction. Also interesting is the low 
ratio for the final states with three kaons, where we do 
not observe a significant signal. These differences in frag- 
mentation could be accounted for by changes in the pa- 
rameters within JETSET, and by the addition of reso- 
nant contributions, but a detailed study of this requires 
a larger data sample. 

We also make a comparison of the ratio of K® to K ~ 
final states in data with the same ratio in Monte Carlo. 
In the range M(X S ) = 1.1 — 2.8 GeV this gives a double 
ratio of 1.00 ±0.21, which is consistent with the assump- 
tion of isospin symmetry. 

The ratios in Table [H] are applied as weights to the 
generated and reconstructed signal Monte Carlo events 
to correct for the observed fragmentation. The reduction 
in the Ktt final states and the increase in the high multi- 
plicity final states reduces the signal efficiencies by 10 % 
to 25%, depending on the M(X S ) bin. The weights also 
lead to an increase in the cross-feed background. 

The selected 38 X s final states do not account for all 
the states produced in B — » X s -f. To obtain the total 
B — > X s j branching fraction in each M (X s ) bin we need 
to correct for the fraction of missing final states, which 
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TABLE II: Ratios of data to Monte Carlo yields for various 
categories of final states. These are used to adjust the frag- 
mentation in the signal Monte Carlo. 



100 



Final States 



Data/Monte Carlo 



is- „o is" 

K TT , K q TT 



0.50 ±0.07 
0.19 ±0.12 



R--K+ 7T _ , K° s 7T+ 7r 
K~TT + 7T°, K° s TV' TT° 



1.02 ±0.14 
1.34 ±0.24 



K TT + TT TT + , Kg TT 
K-TT+ TT- TT°, K?, 7T + TT 



TT TT 





2.67 ±0.96 
1.29 ±0.61 



K TT TT , K s TT TT 

K TT TT TT , K s TT TT TT 



•±1.33 



K TT + TT TT + TT , K s TT + TT TT + TT 
K~~TT + TT~ TT + TT° , K® TT + TT~ TT~ TT° 



i on+1.55 
!.OZ_ 1 32 



K~n, k° s r,, k-t] tt+ 

K° s j] TT~ , K~r) tt°, K° s rj tt° 

K~ Tj TT + TT~ , Kg T\ TT + TT~ 
K~ri TT+ TT°, K° s Tj TT" TT° 



0.83 



+ 1.00 
0.83 



K~ 


K + K~ 


, K~K + K° S 




K~ 


K+K~ 


tt+, K-K+K% 7T- 


27+ ' 54 


K~ 


K+K- 


TT°, K-K+Kl TT° 





increases from 25 % to 70 % for M{X S ) = 0.6 to 2.8 GeV 
(see Fig. . The 25 % fraction missing at all hadronic 
masses comes from K®. This fraction is equal to the K® 
sample, with an uncertainty determined by our isospin 
asymmetry measurement (see section X). 

An analysis of the final states generated in the signal 
Monte Carlo sample shows that the largest missing con- 
tribution comes from high multiplicity final states with 
a kaon and > 5tt. There are also missing contributions 
of a few percent in the highest mass bins from higher 
multiplicity final states with r\ mesons, three kaons and 
baryons. Smaller contributions come from rare radiative 
meson decays and final state radiation. 

We use the results of the fragmentation study to cor- 
rect for the missing high multiplicity states with rj and 
three kaons with the weights found for the observed final 
states. Since the kaon and 4tt modes are consistent with 
the Monte Carlo expectation within a large statistical 
error, we do not adjust the fractions of kaon and > bu 
modes. For these final states, and for missing fractions 
where there is no information from the reconstructed fi- 
nal states, we assign systematic errors which allow the 
missing fraction to vary by a factor two relative to the 
predicted value. The uncertainty in the missing fractions 
is the dominant systematic error in the high mass bins. 
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FIG. 4: Missing fractions with systematic errors as a function 
of the hadronic mass. 



VIII. SYSTEMATICS 

The following systematic errors are independent of 
M(X S ). There is a 1.1 % uncertainty in our knowledge of 
the number of B mesons in our data sample. There is a 
2.5 % uncertainty in the efficiency of the initial selection 
of the high energy photon, but we apply additional cri- 
teria to isolate the high energy photon, as well as vetoes 
if it forms a ir° or rj meson. The efficiency of these addi- 
tional selections is checked using the K*j region as a con- 
trol sample, and we assign an additional uncertainty of 
1.5 %. The total systematic error independent of M(X S ) 
is 3.1%. 

The systematic errors that depend on M(X S ) are sum- 
marized in Table II I II There are several different cate- 
gories of systematic errors associated with detection effi- 
ciencies, fitting, the modeling of peaking backgrounds, 
fragmentation corrections, and the estimation of the 
missing fractions. 

The detection efficiency errors come from our knowl- 
edge of the differences between Monte Carlo and data 
obtained from control samples. They reflect the accu- 
racy of our modeling of the detector response, and in- 
clude tracking efficiency (2%), particle identification for 
charged kaons (1%), reconstruction (3%), and pho- 
ton detection efficiency (2.5%). These differences give 
not only systematic errors but also shifts in the signal 
detection efficiencies, as mentioned in Section The 
fitting errors come from varying the fixed parameters 
in the data fits. The signal peak position in tues is 
varied by ±0.4MeV, and the width by ±0.1 MeV to 
allow for a possible variation as a function of M(X S ). 
The width is varied in a correlated fashion with the a 
and tail parameters of the Crystal Ball function. These 
changes in the signal shape alter the signal yields by be- 
tween and 5 events in each M(X S ) bin. The peaking 
background shape is varied within a range allowed by 
the cross-feed and BB Monte Carlo samples. This gives 
smaller changes in the signal yield of about one event in 
each bin. The continuum shape and normalization are 
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TABLE III: Contributions to the M(X s )-dependent system- 
atic error on the branching fraction from detection efficiency, 
fitting, peaking background, fragmentation and missing frac- 
tions, are shown in % as a function of hadronic mass. The 
total systematic errors also include a 3.1 % systematic error 
that is independent of M(X a ) and not listed in the Table. 



M(X a ) 


Detector 


Fitting 


Peaking 


Fragmen- 


Missing 


Total 


(GeV) 


Efficiency 




Background 


tation 


Fraction 




0.6-0.7 
0.7-0.8 


5.2 
5.3 


+20.8 

-21.4 
+33.1 
— 41 .0 


8.1 
2.5 




+ 1.8 

- 2.0 
+ 1.8 

— 2.0 


+ 23.1 

-23.6 
+33.7 
—41.4 


0.8-0.9 


5.4 


+ 2.3 
— 2.3 


1.1 




+ 1.8 
— 2.0 


+ 6.6 
— 6.6 


0.9-1.0 


5.3 


+ 2.2 
— 2.2 


0.8 




+ 1.8 
— 2.0 


+ 6.4 
— 6.4 


1.0-1.1 


5.3 


+10.2 
— 10.4 


3.2 


13.7 


+ 2.6 
— 2.2 


+ 18.5 
— 18.5 


1.1-1.2 


6.2 


+11.1 

— 11.2 


1.7 


5.4 


+ 3.4 
— 2.5 


+14.5 
— 14.4 


1.2-1.3 


6.4 


+ 5.5 
— 5.6 


1.8 


4.5 


+ 3.4 
— 2.4 


+10.5 
— 10.3 


1.3-1.4 


6.6 


+ 4.3 
— 4.5 


1.8 


4.5 


+ 3.3 
— 2.4 


+10.0 
— 9.8 


1.4-1.5 


6.7 


+ 4.6 
- 4.7 


2.1 


5.8 


+ 3.9 
— 2.8 


+ 11.1 
-10.8 


1.5-1.6 


6.9 


+ 2.1 
- 2.4 


2.9 


1.9 


+ 4.7 
- 3.0 


+10.5 
- 9.9 


1.6-1.7 


7.0 


+ 3.1 
- 3.5 


4.0 


5.3 


+ 6.1 
- 3.9 


+12.0 

-11.1 


1.7-1.8 


7.1 


+ 4.1 
- 4.2 


2.0 


5.5 


+ 7.5 
- 4.5 


+12.7 

-11.2 


1.8-1.9 


7.2 


+ 3.2 
- 3.4 


2.2 


6.5 


+ 10.6 
- 5.9 


+15.0 
-12.2 


1.9-2.0 


7.1 


+32.6 
-32.4 


9.3 


5.1 


+ 15.0 
- 8.0 


+38.2 
-35.8 


2.0-2.2 


7.2 


+ 11.6 
- 9.4 


5.1 


6.2 


+23.4 
-12.7 


+28.3 
-19.2 


2.2-2.4 


7.5 


+ 10.8 
- 9.1 


8.9 


7.1 


+36.2 
-19.8 


+40.2 
-25.8 


2.4-2.6 


7.5 


+ 7.6 
- 7.6 


11.1 


9.3 


+55.2 
-29.7 


+58.1 
-34.8 


2.6-2.8 


7.9 


+68.3 
-96.1 


65.7 


11.7 


+71.3 
-39.9 


+119.0 
-124.0 


0.6-2.8 


6.1 


+ 3.2 
- 3.7 


1.6 


5.9 


+ 13.8 
- 7.6 


+ 16.7 
-12.2 



varied by the difference between the continuum Monte 
Carlo and the off-peak data, but this gives very small 
changes in the signal yields. 

The peaking background normalization error is treated 
separately. In the region M(X S ) < 1.8 GeV the main 
uncertainty comes from the charmless hadronic decays 
B — > X s tt°. While some of these have been mea- 
sured, others are yet to be observed. Combining the 
effects of the measured and the unmeasured modes we 
assign an uncertainty of 50 % to the normalization of 
this contribution in the low M(X S ) region. In the re- 
gion M(X S ) > 1.8 GeV the background is primarily from 
hadronic b — > c decays, such as B — > D*p~. By taking 
the weighted average of the uncertainties on the branch- 
ing fractions of the components of this background, we 
assign an uncertainty of 15 % to the normalization of this 
contribution in the high M(X S ) region. The cross- feed 
background depends on the modeling of the B — ► X s j 
process. We use the fragmentation weights and the mea- 
sured spectral shape to adjust the signal Monte Carlo to 
match the data, correct our prediction for the peaking 
cross-feed background, and assign an uncertainty to this 
correction. 

As discussed in Section |VlI| we have studied the dif- 



ferences in fragmentation between data and signal Monte 
Carlo, and re- weighted the signal Monte Carlo to correct 
for the differences found. The weighting factors are listed 
by category in Table ITT1 We vary the weights by their sta- 
tistical errors to evaluate the fragmentation systematics. 
For the 77 and K+4tt categories the weights are consistent 
with one (with larger errors), and we restrict the range to 
be between 0.5 and 2.0. Each weight is varied separately, 
and assumed to be uncorrelated with the other weights. 
The fragmentation errors are limited by the statistics of 
the data sample. In the bin M(X S ) = 1.0 - 1.1 GeV 
the efficiency is computed from the average of the reso- 
nant and non-resonant model efficiencies, and we take the 
difference between them as an additional fragmentation 
systematic of 12.8% in this bin. 

The estimation of the missing fractions is the largest 
systematic error. It is determined by varying the missing 
fractions within the ranges shown in Figure We vary 
all the uncertainties in the missing fractions together, ei- 
ther increasing them all, or decreasing them all. When 
we do this we adjust the sum of the reconstructed final 
states so that the total B — ► X s j rate is unchanged, 
i.e. we actually adjust the relative proportions of recon- 
structed and missing final states. 



IX. BRANCHING FRACTION RESULTS 

The branching fractions in each hadronic mass bin are 
obtained using the signal efficiencies shown in Figure ^ 
the signal yields given in Table |U the fraction of recon- 
structed final states (which is 1 minus the fraction of 
missing final states shown in Figure 0J, and the total 
number of BB pairs in the sample. The systematic stud- 
ies which affect each of these quantities were discussed in 
Section 

Table lTvl shows the results for the branching fraction in 
each hadronic mass bin, as well as the result for the whole 
mass range. Table Ivl shows the corresponding branching 
fractions in bins of the photon energy using Equation 
Q to translate between hadronic mass and photon en- 
ergy. We have taken into account the different bin sizes 
in transforming between B(M(X S )) and B{E 1 ). The 
corresponding hadronic mass and photon energy spec- 
tra are shown in Figure where theoretical predictions 
are shown which will be discussed in Section 1X11 The 
hadronic mass resolution of 10 MeV converts into an E 7 
resolution of 1-5 MeV, and the K* peak can be clearly 
seen in both the hadronic mass and photon energy spec- 
tra. 



X. ISOSPIN ASYMMETRY 

We define the isospin asymmetry as the ratio: 
X sSl ) - Y{B- -> X sul ) 



A - T{B ° 



T(B° -> X sdl ) + T(B~ -> X sul ) 



(2) 
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TABLE IV: Branching fractions in bins of hadronic mass with 
statistical and systematic errors. The bottom line shows the 
total branching fraction obtained from the separate fit to the 
data over the full M(X S ) range, and not from the sum of the 
individual bins. 



TABLE V: Branching fractions in bins of photon energy with 
statistical and systematic errors. 



M(X S ) 


(GeV) 


B(M(X S ))/100 MeV 


(io-°) 


0.6 


- 0.7 


0.4 


± 


0.5 


+ 0.1 
- 0.1 


0.7 


- 0.8 


0.3 


± 


0.7 


+ 0.1 
- 0.1 


0.8 


- 0.9 


20.8 


± 


1.2 


+ 1.3 

- 1.3 


0.9 


- 1.0 


19.6 


± 


1.4 


+ 1-2 

- 1.2 


1.0 


- 1.1 


3.3 


± 


1.2 


+ 0.6 

- 0.6 


1.1 


- 1.2 


6.2 


± 


2.4 


+ 0.9 

- 0.9 


1.2 


- 1.3 


18.1 


± 


3.4 


+ 1.9 

- 1.9 


1.3 


- 1.4 


27.6 


± 


4.6 


+ 2.8 
- 2.7 


1.4 


- 1.5 


OO c\ 


_1_ 
± 


o.U 


+ 2.5 
- 2.5 


1 K 


- l.O 


29.8 


± 


6.0 


+ 3.1 
- 3.0 


1.6 


- 1.7 


28.0 


± 


7.2 


+ 3.3 
- 3.1 


1.7 


- 1.8 


26.9 


± 


8.1 


+ 3.4 
- 3.0 


1.8 


- 1.9 


40.6 


± 


9.7 


+ 6.1 
- 5.0 


1.9 


- 2.0 


8.0 


± 


11.7 


+ 3.1 
- 2.9 


2.0 


- 2.2 


21.0 


± 


9.6 


+ 5.9 

- 4.0 


2.2 


- 2.4 


26.1 


± 


12.0 


+ 10.5 
- 6.7 


2.4 


- 2.6 


28.0 


± 


16.0 


+ 16.2 
- 9.7 


2.6 


- 2.8 


-3.7 


± 


18.8 


+ 4.4 
- 4.5 


B (10" 6 ) 


0.6 


- 2.8 


327.0 


± 


18.0 


+55.0 
-40.0 



The Standard Model predicts no isospin symmetry- 
breaking from the dominant penguin diagram for B — ► 
X s "f. Isospin symmetry-breaking effects occur at order 
A/ nib in the heavy quark expansion |38|. due to an- 
nihilation contributions from four-quark operators, the 
chromo-magnetic dipole operator and charm penguins. 
For the exclusive decays B — > K*j, the Standard Model 
predicts a positive value of A - between 5 and 10 % , 
but new physics beyond the Standard Model could en- 
hance the isospin breaking effects. Measurements of the 
B — ► K*j isospin asymmetry from BABAR and BELLE 
are consistent with the predictions of the Standard Model 

HI 

We split the 38 modes into charged and neutral B 
decays, re-fit the data, and calculate the separate effi- 
ciencies and total branching fractions. While the sig- 
nal detection efficiencies are almost a factor of two lower 
for the B~ decays, the backgrounds and missing frac- 
tions are symmetric. Comparing the charged and neu- 
tral branching fraction measurements^using the lifetime 
ratio, t(B-)/t(B°) = 1.086 ±0.017 33], and our recent 
measurement of the production ratio of charged and neu- 
tral B events at the T(4S), B° /B~ = 1.006 ±0.048 Hj, 



> 

CD 



O 
O 



o 



o 

03 



CD 

c 

o 
c 

05 

m 



E 1 (GeV) g(g 7 )/100MeV (10' 



2.593 - 


2.606 


3.3 


± 4.0 


+ 0.8 
- 0.8 


2.579 - 


2.593 


1.9 


± 4.9 


+ 0.6 
- 0.8 


2.563 - 


2.579 


129.2 


± 7.2 


+ 8.1 
- 8.1 


2.545 - 


2.563 


108.9 


± 7.6 


+ 6.6 
- 6.6 


2.525 - 


2.545 


16.7 


± 6.2 


+ 3.2 
- 3.2 


2.503 - 


2.525 


28.6 


±11.1 


+ 4.1 
- 4.1 


2.480 - 


2.503 


76.3 


±14.4 


+ 8.0 
- 7.8 


2.454 - 


2.480 


107.8 


±17.9 


+ 10.8 
-10.6 


2.427 - 


2.454 


82.4 


±18.1 


1 n 9 

- 8.9 


2.397 - 


2.427 


101.6 


±20.3 


+ 10.6 
-10.1 


2.366 - 


2.397 


89.5 


±22.9 


+ 10.7 
- 9.9 


2.333 - 


2.366 


81.3 


±24.6 


+ 10.3 
- 9.1 


2.298 - 


2.333 


115.8 


±27.6 


+ 17.4 
-14.1 


2.261 - 


2.298 


21.8 


±31.6 


+ 8.3 
- 7.8 


2.181 - 


2.261 


52.7 


±24.2 


+ 14.9 
-10.1 


2.094 - 


2.181 


60.0 


±27.6 


+24.1 
-15.5 


1.999 - 


2.094 


59.0 


±33.8 


+34.3 
-20.5 


1.897 - 


1.999 


-7.1 


±36.7 


+ 8.5 
- 8.8 
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FIG. 5: The hadronic mass spectrum (a), and the photon en- 
ergy spectrum (b). The data points are compared to theoret- 
ical predictions (histograms) obtained using the shape func- 
tion (solid line) and kinetic (dashed line) schemes. 
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gives the isospin asymmetry over the range M(X S ) — 
0.6 - 2.8 GeV: 

A _ = -0.006 ± 0.058 ± 0.009 ± 0.024 

The errors are statistical, systematic and due to the 
production ratio B Q / B~ , respectively. Most of the sys- 
tematic errors on the branching fractions cancel in the 
ratio. The residual systematic errors that are relevant 
to the isospin asymmetry are ±0.001 from the detec- 
tion efficiency corrections, and +0.004 from a study of 
the fragmentation differences between charged and neu- 
tral B modes in which we allow for different B~ and B° 
weights. The largest contribution to our systematic error 
is ±0.008 due to the uncertainty in the lifetime ratio. 

XI. FITS TO THE SPECTRUM, EXTRACTION 
OF m b AND pi, AND INCLUSIVE BRANCHING 
FRACTION 

In the following section, we evaluate the results in the 
context of recent QCD calculations in the shape func- 
tion 0, El El El El E3 and kinetic |2j schemes. 
From a fit to the spectrum we evaluate the 6-quark 
mass m(,(/i), and the kinetic-energy parameter pl(p>), 
where the reference scales are taken to be p = 1.5 GeV 
in the shape function scheme, and p = 1.0 GeV in the 
kinetic scheme. We have set the chromomagnetic oper- 
ator fiQ = 0.35 GeV 2 , and the Darwin and spin-orbit 
operators to p D = 0.2 GeV 3 and p LS = -0.09 GeV 3 in 
the kinetic scheme |4J|. The photon and hadronic mass 
spectra are equivalent, so we can fit either one to extract 
the heavy-quark parameters. 

We use a x 2 method to fit the spectrum and find the 
best values of the parameters mb and pi, adding an ad- 
ditional constraint on the normalization of the spectrum 
from the value of B(b — ► sj) measured over the full range 
M(X S ) = 0.6 - 2.8 GeV. The fit method takes into ac- 
count the asymmetry of the systematic errors and the 
large bin-to-bin correlations. 

The spectrum is fit using the expected spectral shape 
from the two schemes, except that at low hadronic mass 
we replace the inclusive theoretical distribution with a 
Breit-Wigner to represent the K* resonance. The tran- 
sition point between the K* and inclusive distributions is 
a free parameter of the fit. As a cross-check we have also 
performed fits to the spectrum where we treat the K* re- 
gion M(X S ) = 0.6 — 1.2 GeV as a single bin. In this case 
a Breit-Wigner shape and a transition point are unneces- 
sary to describe the data. The results of this cross-check 
agree with the default fits with the K* included. 

The signal Monte Carlo, which is used to determine 
the signal efficiencies and the cross-feed background, de- 
pends on the mass scheme used and the corresponding 
heavy-quark parameters. Starting from the initial fits 
to the measured spectrum we modify the signal Monte 
Carlo with the fitted parameters, and revise our esti- 
mates of efficiencies and cross-feed. We re fit the data, 
re-calculate the branching fractions and then re-fit the 



spectrum. This procedure leads to small changes in the 
heavy-quark parameters which are well within the errors. 

The final results after the re-fitting are shown in Ta- 
ble ^T] for mt, and /i 2 . In the shape function scheme we 
present results using three different models for the shape 
function E3- The exponential and hyperbolic models 
give very similar results, and the exponential model is 
taken to be the default in this scheme. Slightly differ- 
ent results are obtained with the Gaussian model. In 
the kinetic scheme, we use the two functions provided 
by the authors, and quote the average ^(j- We show 
the fits to the measured hadronic mass and photon en- 
ergy spectrum in Figure [5] from which it can be seen 
that the spectrum is well described and the difference be- 
tween the two schemes is small. The central values and 
error ellipses for the shape function scheme of the fit- 
ted heavy-quark parameters are shown in Figure (the 
points of the exponential model ellipse are given in Ap- 
pendix ^J. For comparison with previous measurements 
we have also fitted the spectrum to the older model of 
Kagan and Neubert g]. We find m b = (4.79t°;^) GeV 
for the b quark mass and A x = (-0.24^°;^) GeV 2 for the 
kinetic parameter, with a 94 % correlation between them. 

TABLE VI: Heavy-quark parameters nib and pi from fits to 
the spectrum using the three different models of the shape 
function scheme and using the kinetic scheme. For each 
scheme the first two columns are from a fit with just statisti- 
cal errors, and the last two columns from a fit with statistical 
and systematic errors. The correlation coefficients between 
the two parameters are -94% and -92% in the two schemes. 



Theoretical Stat. Errors Stat.+Syst. Errors 

Scheme m h (GeV) ^(GeV 2 ) m 6 (GeV) ^(GeV 2 ) 

Shape Function 

Exponential 4.65 ± 0.04 0.19 ± 0.06 4.67 ± 0.07 0.16±g; g 

Hyperbolic 4.64 ± 0.04 0.20 ± 0.06 4.67 ± 0.07 0.17+° Jg 

Gaussian 4.68 ± 0.04 0.12 ± 0.06 4.73 ± £07 °- 07 to'.06 

Kinetic 4.67 ± 0.04 0.32 ± ° ° 7 4 4.70 ± ° ° t 0.29±°;™ 



The inclusive branching fraction is obtained from a 
fit to the mES distribution of data events over the full 
M(X S ) range, corresponding to £L > 1.90 GeV (Fig- 
ure |3J). In the shape function scheme we obtain B(b — > 
s 7 ) = (3.27 ± 0.18t°;!oiao6) x 10 ~ 4 and in the kinetic 
scheme 8(6 -»■ try) = (3.27 ± 0.18±g;lg±g-.;|) x 10" 4 . The 
errors are respectively statistical, systematic and due to 
the variation of the shape parameters. We quote the av- 
erage of the results from the two theoretical schemes: 

8(6 -» s 7 ) = (3.27±0.18l£^tg;^)xl0- 4 , £ 7 > 1.9 GeV 

The branching fraction can be extrapolated to a lower 
photon energy using the fits to the spectrum. Again we 
quote the average of the two schemes: 

B{b -» sj) = (3.35±0.19^i^;^)xl0^ 4 , £ 7 > 1.6 GeV 
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FIG. 6: Error ellipses corresponding to A^ 2 = 1 from the fit 
to the spectrum in the shape function scheme using the expo- 
nential (solid line), hyperbolic (dotted-dashed) and Gaussian 
(dashed line) models of the shape function. Both statistical 
and systematic errors have been taken into account. 



where the small uncertainties from the extrapolations in 
the two schemes enter into the model error through the 
variation of the fitted parameters. 



XII. MOMENTS OF THE PHOTON ENERGY 
SPECTRUM 



The first moment is defined as the average of the pho- 
ton energy spectrum, (i? 7 ), while higher moments are 
defined as {{E 1 — (E 1 )) N ) 1 where N is the order of the 
moment under investigation. The values of the moments 
depend on the range of the photon energy spectrum used 
to calculate them. We vary the range by considering five 
different minimum photon energies E™ m — 1.897, 1.999, 
2.094, 2.181 and 2.261 GeV. These values correspond 
to the boundaries of the highest bins in hadronic mass, 
M(X S ) = 2.8, 2.6, 2.4, 2.2 and 2.0 GeV. 

The results for the first, second and third moments 
of our photon energy spectrum as a function of the 
minimum photon energy are shown in Table IVIII The 
third moments can be used to test the predictions for 
the B — > Xsj decay spectrum by dressed gluon expo- 
nentiation 41], but for now they are statistically lim- 
ited. Figure shows the first and second moments as a 
function of the minimum photon energy, together with 
the predictions from Ref. [2(j, which uses the parame- 
ters m b = 4.61 GeV and n\ = 0.45 GeV 2 that are ob- 
tained from fits to the b — ► civ moments 21] . The 
solid lines represent the band allowed by theoretical un- 
certainties 40] . While the experimental errors decrease 
rapidly as the minimum photon energy is raised, the the- 
oretical errors increase due to the "bias" corrections de- 
scribed in Ref. 20] . The agreement between the b — > s-f 
and b — > civ moments is good, and well within the ex- 
pected theoretical uncertainties. This demonstrates a 
non-trivial consistency between two different classes of 



inclusive b decays. 



TABLE VII: First, second and third moments of the photon 
energy spectrum as a function of the minimum photon energy 
with statistical and systematic errors. 



1,-n't i n 
fc 7 



( GeV) 



{E,) (GeV) 



1.897 
1.999 
2.094 
2.181 
2.261 



2.321 ± 0.038 lo.o38 

+0.014 
-0.029 
+0.007 



2.314 ± 0.023 
2.357 ± 0.017 
2.396 ± 0.013 
2.425 ± 0.009 



+0.003 



+0.002 



(GeV) {(En/ — (E-,)) 2 ) (GeV 2 ) 



E", 



1.897 
1.999 
2.094 
2.181 
2.261 



0.0253 ± 0.0101 laoo28 



0.0273 ± 0.0037 



+0.0015 



0.0183 ± 0.0023 ±°'°ooo° 
0.0115 ± 0.0014 laooos 
0.0075 ± 0.0007 



(GeV) ((E-y — (Ej)) 3 ) (GeV 3 ) 



1.897 
1.999 
2.094 
2.181 
2.261 



-0.0006 ± 0.0012 ±0.0002 

-0.0009 ± 0.0006 1o.oot4 

-0.0005 ± 0.0003 toloooi 

-0.0001 ± 0.0001 toloooo 

+0.0001 ± 0.0001 Innnnn 



We perform fits to the first and second moments to 
obtain values of the heavy-quark parameters which are 
less dependent on the details of the spectral shape. We 
take into account the correlation coefficients between the 
errors on the moments. The full correlation matrices are 
given in Appendix lAl 

Using the shape function scheme, but ignoring theo- 
retical uncertainties, we fit to the first and second mo- 
ments at the lowest minimum photon energy E™ m — 
1.897 GeV. The fitted values are m b = (4.60± o ;^) GeV 
and ^i 2 = (0.19±q Jo) GeV 2 . These are in agreement 
with the fit to the full spectral shape, but have larger 
errors. A fit to the moments at the higher photon en- 
ergy E™ in = 2.094 GeV gives m b = (4.531^;^) GeV and 
[i\ = (0.351^4) GeV 2 , and a fit to the highest photon 
energy E™ n — 2.261 GeV reproduces the results from 
the fit to the spectrum. 

Results with similar precision are obtained from fits 
using the kinetic scheme if theoretical uncertainties are 
ignored. At the higher photon energies the theoretical 
predictions for inclusive quantities such as moments are 
not very reliable, as shown by the solid lines in Figure 
It appears that an optimal combination of experimental 
and theoretical uncertainties occurs in the range EV? m = 
1.9 - 2.1 GeV. 
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E™(GeV) 

FIG. 7: First (a), and second (b) moments as a function of the 
minimum photon energy. The dotted lines show the predicted 
central values based on fits to the 6 — > civ moments pl|. and 
the solid lines the theoretical uncertainties from the kinetic 
scheme HJ|4(J. 



XIII. CONCLUSIONS 

We have measured the inclusive b — > sj branch- 
ing fraction for E 7 > 1.9 GeV and extrapolated it to 
Ej > 1.6 GeV. Our result is in good agreement with the 
world average. Although our measurement is currently 
systematics limited, with more data we expect to improve 
our understanding of the missing fraction of high multi- 
plicity X s hadronic final states, and therefore reduce our 
dominant systematic error. We have made the first mea- 
surement of the isospin asymmetry between B~ — » A sfi 7 
and £>° — > X s( ff, and found that it is consistent with zero 
within the experimental uncertainty, which is mainly sta- 
tistical. 

We have made a measurement of the b — > spectral 
shape over the range E 1 > 1.9 GeV. After taking into 
account the presence of the if* (892) resonance, the shape 
is found to agree well with two recent theoretical calcula- 
tions. Fits to the spectrum are used to give values for the 
heavy-quark parameters m& and /i^ in the two schemes. 
We calculate the first, second and third moments of the 
photon spectrum for five different minimum values of the 
photon energy between 1.90 and 2.26 GeV. The moments 
are in good agreement with predictions based on fits to 
the measured b — > civ moments. Fits to the photon en- 
ergy moments give heavy-quark parameters which agree 
with the fits to the full spectrum, but are less accurate. 
Information from b — > 57 and b — > civ can be combined 



to obtain tighter constraints on mb and /ij . This will lead 
to improved extractions of \V u b\ from the measurements 
of b — > ulv decays. 

While this paper was in preparation, new references 
appeared concerning both the photon energy spectrum 
and the moments. Ref. |4lj computes the photon energy 
spectrum and moments by resummed perturbation 
theory, using the technique of Dressed Gluon Exponen- 
tiation. The predicted spectrum extends smoothly into 
the non-perturbative region and tends to zero at the 
physical endpoint. The predictions for the moments have 
been compared to our results and are consistent with 
them. Ref. 0] presents fits to the BELLE spectrum |24| 
in the shape function scheme. Ref. 0] shows predictions 
for the photon energy moments using an OPE approach 
calculated to NNLO accuracy. This does not require 
spectral information from the shape functions. These 
predictions are fit to the experimental moments from 
the BELLE spectrum and to the moments from this 
analysis. The new preprints indicate the current interest 
in the extraction of heavy-quark parameters from the 
shape of the 6 -> S7 spectrum, and show how our 
data will contribute to improved knowledge of these 
parameters. 
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APPENDIX A 

In this appendix we present more detailed information 
concerning the fit to the spectrum and the correlation 
matrices for the moments. 



Ellipse from the Fit to the Spectrum for the 
Shape Function Scheme 



TABLE IX: Statistical correlation coefficients for the mo- 
ments with different minimum cuts on the photon energy. 



Table IVTlTl lists the points on the ellipse from the fit to 
the spectrum using the exponential model in the shape 
function scheme and shown in Figure 



TABLE VIII: Heavy-quark parameters, mb and fj%, evalu- 
ated along the A\ 2 = 1 contour resulting from the fit to the 
spectrum using the shape function scheme. 



mi,(GeV) /J, 


^(GeV) 


m b (GeV) 




4.618 


0.249 


4.730 


0.088 


4.629 


0.240 


4.721 


0.095 


4.636 


0.233 


4.716 


0.099 


4.642 


0.227 


4.709 


0.106 


4.647 


0.222 


4.702 


0.112 








nils 


4.662 


0.205 


4.692 


0.124 


4.669 


0.197 


4.686 


0.129 


4.675 


0.190 


4.681 


0.135 


4.680 


0.183 


4.677 


0.140 


4.685 


0.177 


4.672 


0.146 


4.690 


0.171 


4.667 


0.152 


4.694 


0.165 


4.659 


0.162 


4.702 


0.154 


4.650 


0.173 


4.709 


0.145 


4.642 


0.183 


4.715 


0.136 


4.634 


0.194 


4.721 


0.127 


4.627 


0.204 


4.726 


0.120 


4.620 


0.214 


4.730 


0.113 


4.614 


0.224 


4.734 


0.106 


4.608 


0.233 


4.737 


0.100 


4.604 


0.241 


4.739 


0.095 


4.601 


0.249 


4.741 


0.090 


4.601 


0.255 


4.742 


0.087 


4.603 


0.257 


4.740 


0.084 


4.607 


0.257 



E mm (GeV) 1.897 1.999 2.094 2.181 2.261 


(E) 




1.897 


1 0.55 0.28 0.12 0.033 




1.999 


1 0.47 0.19 0.046 


(E) 


2.094 


1 0.48 0.18 




2.181 


1 0.45 




2.261 


1 


({E-(E)f) 




1.897 


1 0.28 0.03 -0.07 -0.06 




1.999 


1 0.16 -0.14 -0.14 


((E-{E)f) 


2.094 


1 0.13 -0.09 




2.181 


1 0.17 




2.261 


1 



(E) 



1.897 
1.999 
2.094 
2.181 
2.261 



((E-(E)f) 

-0.90 -0.41 -0.23 -0.12 
-0.79 -0.40 -0.19 
-0.79 -0.42 
-0.82 



-0.05 
-0.08 
-0.18 
-0.37 
-0.76 



TABLE X: Systematic correlation coefficients for the mo- 
ments with different minimum cuts on the photon energy. 



E mm (GeV) 1.897 1.999 2.094 2.181 2.261 

(E) 



2. Correlation Matrices for the Moments 

Here we present the correlation matrices for the first 
and second moments. Table llXl shows the statistical cor- 
relation matrix and Table Ixl shows the systematic corre- 
lation matrix. 





1.897 


1 0.98 


0.97 0.93 


0.85 




1.999 
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0.98 0.95 


0.84 


(E) 


2.094 




1 0.97 


0.86 




2.181 
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0.92 




2.261 






1 






{(E 


- (E)) 2 ) 






1.897 


1 0.87 


0.86 0.62 


0.30 




1.999 
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0.90 0.70 


0.27 


((E-{E)f) 


2.094 




1 0.85 


0.38 




2.181 




1 


0.67 




2.261 






1 



((E-(E)f) 



(E) 



1.897 
1.999 
2.094 
2.181 
2.261 



-0.93 -0.87 -0.91 -0.71 
-0.90 -0.91 -0.74 
-0.90 -0.75 
-0.73 



-0.26 
-0.22 
-0.21 
-0.20 
-0.08 



